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In this letter we investigate the anisotropic g-state Potts model on the square 
lattice]!], |2|. This model is defined by the Hamiltonian 

H — — ^2(J x S(o- Xj y - cr x+lt y) + J y S(a Xt y - CTjj^+l)) (1) 

where the variables a X}V may take one of q distinct values, and (x, y) runs over the 
coordinates of the sites of the lattice. 

Much is known in the cases when J x and J y are of the same sign0, [3], corresponding 
to the ferromagnetic Potts model (J x > 0, J y > 0) and the antiferromagnetic Potts 
model (J x < 0, J y < 0)(4|. Much less is known in the mixed interaction case, J x > 
and J y < 0. This will be the object of this letter. 

There have been a number of attempts to study the mixed interaction model, 
particularly for the case q = 3[]5j, |6|, [7], Kinzel, Selke and Wu0 obtained a candidate 
critical line by applying a one step Migdal-Kadanoff procedure, thereby mapping the 
model onto the ferromagnetic Potts model for which the exact expression for the critical 
line is known. This gave the critical line as 

(1 + exp((3J x ))(l - exp(pjy)) = q. (2) 

Spurred on by the fact that this line is exact for the exactly known q = 2 case, it was 
suggested that perhaps this equation may remain exact for other values of q. They also 
gave numerical results from Monte Carlo simulations which suggest a transition line at 
lower temperatures than given by @. The transition line was estimated, however, by 
extrapolating the inverse correlation lengths by hand to zero for lattice sizes less than 
about 60. 

The possibility that (0) may correspond to an exact result was refuted by the 
observation that, except for q = 2, this line does not remain invariant under the 
symmetries of the model0. The symmetries of the model may be used to construct 
a symmetry group 0. The critical line must remain invariant under the action of 
these symmetries. Using the group of symmetries of the model, Truong|| constructed 
extended duality relations which he used to look for candidate critical lines. For the 
3-state Potts model this gave: 

exp(P(J x + J y )) + 2 eM(3J x ) - exp(pjy) + 1 = 0. (3) 

Interestingly, this line has the same functional form as the ferromagnetic transition line 
but with J y — > — J y . This line is indeed invariant under all the symmetries of the model, 
but does not seem to correspond to the transition line when compared with numerical 
results of Kinzel et al 0, Yasumura|| and of the present work. This leaves the 
interesting question: Does this very special line have a corresponding physical meaning? 
This question remains open. 

It has been conjectured that the transition for q = 3 is of Kosterlitz-Thouless type. 
This conjecture is based on the free-fermion approximation applied to a related clock- 



model[10], and supported by numerical results on a related one dimensional quantum 



Potts model||TI|. In this letter we present convincing evidence that the transition for 
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q = 3 is indeed of the Kosterlitz-Thouless type. The critical exponent v is given as a 
function of q and shown to diverge at q — 3. The critical temperatures are also given, as 
a function of q, for the case J x = —J y > 0. No transition is found for q > 3. This work 
provides a more accurate estimate of the critical line for q = 3, which is substantially 
lower than previous estimates |5], [7|, |8|J. 

The numerical results were obtained using the numerical Transfer Matrix Method 
coupled with phenomenological renormalisation||i"2| . The method consists of expressing 
the partition function of an infinite strip of finite width L in terms of a matrix product. 
Imposing periodic boundary conditions in the transfer direction (along the direction 
taken to be infinite), the partition function for the strip may be written 

Z L = lim TrT^. (4) 

iV-KX) 

The thermodynamic limit corresponds to the limit L —> oo. For integer values of q > 1 it 
is possible to write the transfer matrices directly in terms of spins. Full details are given 
in Foster, Gerard and Puha|T3|]. For non-integer values of q it is necessary to use the 
Kasteleyn-Fortuin mapping[i"2|, which gives the Potts model in terms of bond-clusters 
on the lattice. Full details on the construction of the transfer matrices for this case are 
given in Blote and Nightingale . 

It may be shown that the correlation length in the transfer direction is given by 

& = (5) 

Assuming scale invariance at the critical point, finite size estimates of the critical 
temperature may be identified with solutions of the equation 



L U 

The critical exponent v may be estimated at the fixed points, solutions of (§), by 

i jgg(t/g£ 



(6) 



vljj 



-1 (7) 



Applying these results to the case J x = —J y > gives the results shown in figure [3] 
for the critical temperature as a function of q. In figure |l] (A) the transfer direction is 
taken in the x-direction, i.e. parallel to the ferromagnetic interactions, and in figure [3] (B) 
the transfer direction is taken in the ^-direction, i.e. parallel to the antiferromagnetic 
interactions. In the first case the largest value of q for which a solution was found, 
Qmax(L), decreases with L. On the other hand, g max (L) increases with L in the second. 
This indicates a value of q, g max , such that there exists a critical phase transition for 
(1 < <?max, but not for q > g max - Figure ^] shows the values of g max found in the two 
cases. The curves both extrapolate plausibly to g max = 3 in the thermodynamic limit. 
Using the upper line in figure ^ to extrapolate, we qive g max = 3.00 ± 0.003. This is in 
contradiction with Yasumura||, who gives non-zero values of T c even for q > 3. 

The finite-size estimates of v are shown in figure |] for a transfer direction in the 
x-direction. They converge well to v = 1 for q = 2, and diverge for q = 3. Table [I] 
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Figure 1. The critical temperature as a function of q for a transfer direction along 
(A) the ^-direction and (B) the y-direction. Periodic boundary conditions are taken 
in both directions. 
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Figure 2. The values of q max , the largest values of q for which an estimate of a 
critical temperature could be found, as a function of 1/L. The upper line corresponds 
to the transfer direction in the x-direction, while the lower line of points corresponds 
to the transfer direction in the y-direction. The point q — 3 and 1/L = is shown for 
reference. 



^max 3.5 




shows the finite-size estimates of T c for selected values of q. The estimates of T c are 
extrapolated using a Burlich and Stoer algorithm for q = 1, 1.5,2 and 2.8, but had to 
be extrapolated graphically for q = 3. 

The phase diagram for q = 3 is shown in figure |j. The critical line estimates 
obtained are compared with the exact line for J y /J x > 0, and shown to converge well, 
and are compared to both proposed critical lines. The convergence is much slower for 
Jy/Jx < 0, but clearly rules out both proposed analytic results^, [| and the numerical 
results||, g |. 

In summary, in this letter we have addressed numerically the nature of the criticality 
in the only sector of the anisotropic Potts model for which exact solutions do not as 
yet exist: the ferro/antiferromagnetic Potts model. We have given critical temperature 
estimates and estimates of v as a function q for J y /J x = —1. We find convincing 
evidence that there exists a g max around q = 3 above which there is no transition. At 
Qmax the value of v is diverging, consistent with a Kosterlitz-Thouless type transition. 
This type of transition has already been conjectured for q = 3 using the free-fermion 
approximation |10[| . We conjecture, based on this and our numerical data, that g max = 3. 
This picture is analogous with what happens in the 0(n) spin models in two dimensions, 
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Figure 3. Finite size estimates of v as a function of q. 




Figure 4. Phase diagram for q = 3. The cross shows the estimated critical 
temperature of Kinzcl et al for J y /J x = — 1, along with error bar. The star shows 
the extrapolation of our estimated critical temperatures, also for J y /J x = — 1, with 
error bar. 
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Table 1. Estimates of T c for selected values of q. The transfer direction is taken in 
the x direction. 



L/V 


q = i 


q = 1.5 


q = 2 


q = 2.8 


g = 3 


2/4 


2.701840 


1.677124 


1.172009 


0.742794 


0.670151 


4/6 


2.728013 


1.678307 


1.143901 


0.642965 


0.536430 


6/8 


2.737827 


1.680590 


1.137474 


0.606052 


0.472832 


8/10 


2.742918 


1.682553 


1.135844 


0.590809 


0.438814 


10/12 


2.745847 


1.683806 


1.135255 


0.583354 


0.416617 


CO 


2.75 ±0.01 


1.69 ±0.01 


1.13 ±0.01 


0.57 ±0.01 


0.32 ±0.03 



where there is a transition for — 2 < n < 2, the transition at n = 2 corresponding to the 
Kosterlitz-Thouless transition |T7] . 
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